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ABSTRACT 

Computational methods have been developed and 
successfully used for determining the optimum 
distribution of space radiation shielding on geo- 
metrically complex space vehicles. These methods 
have been incorporated in computer program SWORD 
which uses the full capability of state-of-the-art 
methods for dose evaluation in complex geometry, 
and iteratively calculates the optimum distribution 
of (minimum) shield mass satisfying multiple acute 
and protracted dose constraints associated with 
each of several body organs. The unique and 
effective technique used to accommodate multiple 
constraints, eliminates the awkward discontinuities 
associated with the formulation of inequality con- 
straints, and produces a result meeting mathemat- 
ical tests for optimality. 


INTRODUCTION 


The capability to compute space radiation 
doses at specified points within a vehicle and 
shield of fixed geometry is provided by a number 
of computer programs, e.g. , SIGMA (Ref. 1), 

MEVDP (Ref. 2), and LSVDC4 (Ref. 3). An additional 
capability is provided by techniques developed at 
the McDonnell Douglas Astronautics Company and 
incorporated in the SWORD program (Ref. 4), this 
being the automated computation of the optimal 
shield mass distribution that meets a set of 
radiobiological dose criteria associated with the 
specified vehicle configuration and mission pro- 
file. It is emphasized that this optimization 
function is accomplished without necessitating any 
simplifying assumptions regarding geometrical 
framework, r ad la tion transmission evaluation, etc. , 
that are not also commonly invoked in programs 
which perform dose evaluation only. In fact, SWORD 
incorporates the dose analysis framework of the 
SIGMA program referred to above. To a great extent, 
the capacity of SWORD to accommodate a variety of 
complicating factors influencing shield mass dis- 
tribution is due to the efficacy, yet simplicity, 
of the optimization technique employed. 

The role of the SWORD program in performing 
space radiation shielding analyses is indicated in 
Figure 1, together with those of other space 
radiation analysis programs also developed at the 
McDonnell Douglas Astronautics Company. Overall 
capabilities are summarized in Figure 2. SWORD 
uses basic dose transmission data, in the form of 
one-dimensional point kernel functions, typically 
calculated by the CHARGE program (Ref. 5) for the 
total space radiation environment defined by the 
OGRE program (Ref. 6) for the specified mission. 
These dose transmission data are applied in con- 
junction with ray tracing computations performed 
on a generalized quadric surface representation 


of the vehicle, to compute dose levels to specified 
items, usually the critical organs of crew members. 
The derivative of total dose from all radiation 
sources, with respect to the thicknesses of speci- 
fied candidate shield regions located on various 
surfaces of the vehicle, is also computed. This 
dose derivative information is processed by SWORD 
in an iterative procedure to determine the optimal 
distribution of shield material among such candi- 
date locations as wall structure, biowell for 
solar cosmic ray protection, and personal shields. 
The optimization technique employed is based on a 
particular formulation of the Lagrange multiplier 
constraint equations. Shield shaping over extended 
surface areas is accomplished by subdividing them 
into a number of smaller areas over which shield 
thickness is uniform. SWORD can treat the effect 
of (1) multiple dose constraints (separate con- 
straints for each organ), (2) time- dependent 
astronaut locations (the work-rest cycle influence), 
(3) organ-dependent quality factors or dose distri- 
bution factors, and (4) direct and scattered neutron 
and gamma radiations from on-board nuclear power 
sources. The geometric framework, numerical 
integration schemes, and optimization procedures 
are sufficiently flexible and efficient to allow 
the analysis of any practical space vehicle 
configuration . 


SPACE RADIATION SHIELDING PROGRAMS 


19533 



TRAJECTORY 
MISSION DURATION 
LAUNCH DATE 


OGRE DATA 
SHIELD MATERIALS 


SYSTEM GEOMETRY 
MAN-MODEL GEOMETRY 
WORK-REST CYCLE 
CHARGE DATA 


VEHICLE GEOMETRY 
MAN-MODEL GEOMETRY 
NUCLEAR SOURCE GEOMETRY 
WORK- REST CYCLE 
CHARGE DATA 

CANDIDATE SHIELD LOCATIONS 
DOSE DESIGN CRITERIA 
PHYSICAL SHIELD CONSTRAINTS 


MISSION- INTEGRATED RADIATION 
ENVIRONMENT DEFINED IN INTENSITY, 
ENERGY AND TIME 


DOSE TRANSMI SSION THROUGH 
SPHERICAL SHIELDS 


TIME- INTEGRATED DOSE FOR REALISTIC 
SYSTEM GEOMETRY; SEVERAL PARAMETRIC 
OPTIONS 


TIME- INTEGRATED DOSE FOR REALISTIC 
SYSTEM GEOMETRY (SPACE AND ONBOARD 
NUCLEAR SOURCES). OPTIMIZED SHIELD 
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FIGURE 1 


* The computational techniques described herein were developed by the McDonnell Douglas Astronautics 
Company under the Independent Research and Development Account No. S. 0. 80205-007. 
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SUMMARY OF SWORD CAPABILITIES 


• DESCRIBES SPACE VEHICLE IN A GENERALIZED QUADRATIC SURFACE GEOMETRY 


• USES TABULAR DOSE-SHIELD THICKNESS DATA FROM DETAILED 
SHIELDING PHYSICS CALCULATIONS 

• EVALUATES Ml SSI ON- INTEGRATED CRITICAL ORGAN DOSE USING FLEX I BLE 
NUMERICAL ANGULAR INTEGRATION SCHEME 

• CALCULATES OPTIMAL DISTRIBUTION OF SHIELD MATERIAL SUBJECT TO MULTIPLE 
DOSE CRITERIA 

• USES EFFICIENT, ACCURATE OPTIMIZATION TECHNIQUE 


CRITICAL ORGAN DOSE EVALUATION 

• GEOMETRY MODEL-GENERALIZED QUADRIC SURFACES 

G(X, Y, Z) - Ao+AjK+AgY^ A 3 Z+A 4 X 2 +A 5 Y 2 +A 6 Z 2 + A 7 XY+A 8 YZ+A 9 ZX - 0 

• SPACE RADIATION DOSE-EVALUATED FOR EACH SOURCE AND SUMMED OVER 
VARIOUS TIMEWEIGHTED ASTRONAUT LOCATIONS 


5- +--i 

IQ 

* Klpl IS TABULAR, INTERPOLATED BY POWER LAW BETWEEN ENTRIES 

4tt\ 3/\ 3 / , e Mhi d Tj 

• DATA ARE SAVED FOR EACH RAY DEFINING ITS SOLID-ANGLE WEIGHT, MASS 
TRAVERSED^ . SHIELDS CROSSED, AND COSINE OF ANGLE AT EACH CROSSING 


• DISTRIBUTES SHIELD MATERIAL AMONG NESilD SURFACES OF VEHICLE, BIOWELL, AND 
PERSONAL AND PORTABLE SHIELDS 

. SHAPES SHIELD THICKNESS OVER EXTENDED SHIELD SURFACE AREAS 


> NUCLEAR RADIATION DOSE-EVALUATED CONCURRENTLY WITH SPACE RADIATION 
DOSE EVALUATION USING POINT KERNEL METHOD CALCULATING SCATTERED DOSE 
FOR EACH RAY TRACED IN 4n INTEGRATION AS WtLL AS DIRECT DOSE 


FIGURE 2 


FIGURE 3 


DOSE AND SHIELD MASS EVALUATION 

The significance of an optimum solution to a 
given space radiation protection problem is 
dependent on a number of factors. These include: 

(1) the accuracy with which dose and shield mass 
data, and their derivatives with respect to the 
variables of the problems, can be determined, 

(2) the suitability of the variables for character- 
izing the shield mass distribution, and (3) the 
extent to which all relevant constraints are 
suitably imposed and satisfied. 


Dose Evaluation 

The procedures for dose evaluation that have 
been incorporated in SWORD are outlined in Figure 
3. They have been adapted from the SIGMA code and 
are representative of the state-of-the-art of space- 
radiation-dose analysis. The geometric description 
of the vehicle structure, equipment, stores, fixed 
shields, astronauts, etc., is accomplished using 
generalized quadric surfaces which form the bound- 
aries of homogeneous, contiguous, non- over lapping 
regions. Doses are evaluated at any n\*ber of 
points representing the locations of radiation- 
sensitive body organs by repetitively tracing rays 
from each dose point through the surrounding 
materials, employing the resultant data on mass 
distributions along each ray to determine a 
differential dose contribution for the path. 

Total, mission-integrated critical organ doses 
are obtained by integrating the differential dose 
contributions over solid angle (d0 dy) about each 
dose point and over various time-weighted astro- 
naut locations within the vehicle. The integration 
is accomplished numerically using Simpson's Rule; 

Ai and ^ in Figure 3 are the weighting factors 
associated with the particular values of the 
variables (e, u) • 


The specific dose transmission kernels used 
in an analysis are conveniently supplied in tabular 
form for an arbitrary number of distinct, mission- 
dependent sources (e.g., trapped protons, trapped 
electrons, bremsst rah lung, solar cosmic rays); the 
data are Interpolated by power law between entries 
to obtain dose levels for various values of equiv- 
alent shield thickness. For virtually all analyses 
it is satisfactory to account generally for the 
material-dependence of the dose attenuation kernels 
by applying mass density scaling factors based on 
known radiation-interaction properties. This 
approach is satisfactory for secondary nucleon 
dose as well as for primary proton dose (Ref. 6). 
Since the origin of the dose attenuation kernels 
is Irrelevant to the operation of SWORD, it is 
not restricted to use with a specific one-dimen- 
sional, dose-analysis computer program. 

For each ray traced in the initial evaluation 
of mission-integrated critical organ dose, data 
are stored to facilitate efficient re-evaluation 
of the dose associated with the ray as shield 
thicknesses are iteratively modified. These data 
Include the equivalent thickness of all fixed 
regions traversed, the solid angle weight of the 
ray (incorporating Simpson's rule coefficients 
for the particular values of polar and aslmuthal 
angles defining the ray) , the Indices of all 
shields crossed, and the cosine of the angle at 
the crossing. Information is also saved regarding 
the neutron and gamma differential doses associated 
with the ray; such data represent the dose either 
transmitted directly from an on-board nuclear 
power source or single-scattered in regions lying 
along the ray. The scattered dose calculation is 
limited to scattering in regions external to the 
outermost shielded volume, based on the assumption 
that, for scattering in materials within inhabited 
regions to be significant, the unscattered dose 
would be prohibitive. 
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Shield System Variables 


A framework defining candidate shield 
locations is superimposed on the basic vehicle 
geometry by designating that shield material be 
placed, as required, at a specific bounding sur- 
face of a specific region. Shield material is 
assumed to be uniformly distributed over such 
shield areas, such that a single value of the 
shield thickness measured along a surface normal, 
characterizes the mass of each such candidate 
shield. The relationship of shield mass and the 
derivative of shield mass with respect to the 
thicknesses are then quite simple: 

„ S A t • IS - A 
^ A i i ' 3t t i 


In SWORD the areas of each shield, Ai, may be 
specified by input or may be estimated from data 
obtained during ray tracing computations associ- 
ated with dose evaluation. The latter is con- 
venient for complex shield shapes, but the former 
is usually more accurate. 

This technique for defining shield variables 
is summarized in Figure 4 and illustrated in Figure 
5. In the illustration, a total of seven variables 
have been used to define a shield system consisting 
of an external shell completely enclosing inhabited 
areas (t 1# t 2 , and an internal biowell for 
solar flare protection (t 4 , t 5 , t^, t ? ). Several 
more shields could have been specif lea if, for 
example, it had appeared potentially rewarding to 
use non-uniform shield thicknesses over large areas 
such as those spanned by shields 1, 2 and 3. In 
the actual problem from which this illustration 
was taken, however, the location of fixed equip- 
ment along the vehicle walls was not well defined, 
and it was necessary to assume it to be uniformly 
distributed; in such a circumstance, total shield 
mass is best estimated by also assuming a uniform 
distribution of shield mass over the same areas. 
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VARIABLES AND CONSTRAINTS 

VARIABLES 

• CHARACTERISTIC THICKNESS (g/cm 2 ) AT CANDIDATE SHIELD LOCATIONS 

• THICKNESSES MEASURED ALONG SURFACE NORMALS 
•SHIELD WEIGHT - THICKNESS X AREA 

• TRIPLY-NESTED SHIELDS ALLOWED 


CONSTRAINTS 


• MAXIMUM ALLOWABLE DOSE TO CRITICAL BODY ORGANS (I.E., BFO, SKIN, EYES) 

DOSE FOR PROTRACTED EXPOSURE COMPARED WITH Ml SSI ON- INTEGRATED 
DOSE (ORGAN DOSE SUMMED OVER VARIOUS TIME-WEIGHTED ASTRONAUT 
LOCATIONS) 

DOSE FOR ACUTE EXPOSURE COMPARED WITH DOSE FROM A MAJOR 
FLARE (ASTRONAUT LOCATION MAY BE RESTRICTED DURING EVENT I. E. 

TO A BIOWELL) 

• MINIMUM SHIELD THICKNESS VALUES (>0) 

• MAXIMUM SHIELD THICKNESS VALUES CORRESPONDING TO VOLUME OR WEIGHT 
RESTRICTIONS AT ANY CANDIDATE SHIELD LOCATION 


FIGURE 4 


VEHICLE CALCULATIONAL MODEL 
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O CYLINDER A CONE ® SPHERE □ PLANE 
SH I ELDS DESIGNATED BY Tj 


FIGURES 


In the specification of shield variables for 
a given problem, two capabilities not illustrated 
by Figure 5 are available: (1) shields can be 

tr i P /o7 ne8ted ( doubl y~ ne *ted being illustrated), 
and (2) shields can be portable. Both capabilities 
might be employed if, for example, an astronaut 
wore a personal portable shield for portions of 
the mission, including some time in the biowell. 
Then the effectiveness of such a shield would be 
incorporated in dose and dose derivative calcula- 
tions at all work and rest stations by describing 
the shield geometry at each man-model location, 
but counting its weight only once. 


Constraints on Mass Distribution 

Two types of constraints on shield mass 
distribution can be explicitly included: those 

dealing with maximum levels of radiation exposure 
and those dealing with minimum and maximum values 
of the variables (shield thicknesses). Several 
other constraints are potentially significant, 
but are of such a nature that they can be expressed 
within the overall framework; i.e., restriction to 
non-zero shield thicknesses, restriction of crew 
movement during solar flares, anisotropic radiation 
leakage from on-board nuclear power sources. 

The constraints which can be Imposed are 
listed in Figure 4. Those constraints pertaining 
to radiation exposure are expressed by allowing 
both acute and protracted dose criteria for each 
of several radiation sensitive body organs. These 
dose criteria are understood to be mission-inte- 
grated criteria and are compared, at each step of 
an optimization computation, with dose values 
determined by summing dose contributions for all 
crew stations, explicitly incorporating time 
factors expressing the exposure to each radiation 
source at each station. 
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OPTIMIZATION METHOD 


Numerous mathematical techniques exist for 
computing optimal solutions to engineering prob- 
lems. For many problems, the solutions can be 
obtained directly by solving the set of simulta- 
neous equations for the values of the variables 
satisfying the constraints imposed. In other 
cases, and particularly for geometrically 
realistic shielding problems, the complexity of 
the functional relationships is such that the 
equations can only be solved by iteration on the 
many variables, hopefully converging on a solution 
that is optimum and that can be ascertained for a 
realistic expenditure of engineering labor and 
computer time. 


Capabilities Required 

Among the approaches to optimization of 
radiation shield systems, per se, those that have 
been reported have dealt predominantly with 
optimization of shields for nuclear power systems, 
rather than with shields for space radiations. An 
application of the gradient nonlinear programming 
technique to optimization of a divided shield 
system for a nuclear powered aircraft was reported 
in References 8, 9 and 10; the problem was consid- 
erably simplified by the assumption that the 
source could be adequately represented by a point 
source. Another approach to the same optimization 
problem was reported in Reference 11, in which 
case a restricted representation of shield system 
variables was used to facilitate Iteration on the 
shield thicknesses and on the two Lagrange multi- 
pliers associated with the neutron and gamma dose 
constraints. A shield synthesis technique was 
reported for optimization of compact power reactor 
shields in Reference 12 and an application of the 
technique to optimization of proton shields was 
presented in Reference 13; both applications 
restricted shield geometry to convex shapes. 

Such shield optimization approaches generally 
relaxed one or more important aspects in order to 
facilitate determination of a minimum weight 
solution, i.e. , the detail with which the optimi- 
zation problem can be characterized, the number 
and nature of the constraints which can be simul- 
taneously imposed, or the accuracy of dose and 
dose derivative evaluations. The usual difficulty 
was that the optimization formulation did not have 
the capacity to allcw for much detail in either 
problem geometry (basic system and shield) or 
design criteria. In developing SWORD, however, a 
determined attempt was made to avoid such simpli- 
fications; rather, a number of capabilities 
required of the overall program were identified 
as requirements to be satisfied by the optimization 
procedure. These Included: 

o Compatibility with a detailed geometric 
description of the vehicle, equipment, 
stores, etc., preferably the same 
generalized geometry framework now 
commonly employed in space radiation 
dose evaluation programs and with dose 
transmission data provided by detailed 
physics computations. 


o A flexible system for defining shield 
system variables in a framework which 
can be superimposed on the detailed 
fixed geometry of the vehicle and its 
contents . 

o Capability to accommodate multiple 
radiation level criteria with no a 
priori knowledge of their inter- 
dependence or independence. The 
latter implies that the design values 
for some constraints may necessarily be 
exceeded in meeting the specified values 
for others. Also the relative importance 
of each constraint should be ascertained 
automatically and its effect diminished 
automatically if it becomes an inequality 
constraint. 

o Internal calculation of radiation levels, 
shield weights , and their derivatives , to 
minimize data handling and to permit a 
detailed re-evaluation of all dose and 
weight values as required throughout the 
iterative operations. 

o Procedures for controlling values of the 
variables during iterations so as to 
minimize the quantity of data to be up- 
dated at each step, as this is potentially 
time consuming. 


Candidate Techniques 

Two optimization techniques potentially able 
to satisfy these requirements were identified. 

One consisted of a conventional formulation of a 
multiple-constraint problem using the Lagrange 
multiplier technique and used an iterative proce- 
dure suggested by Arrow and Hurwicz (Ref. 14); 
this method is outlined in Figure 6. The other 
method, developed and applied at McDonnell Douglas, 
reduced the multiple-constraint problem so that it 
can be handled operationally as a single- constraint 
problem. It used a single approximate, continuous, 
functional combination of the multiple constraints, 
which forces the result to converge to the solution 
for the exact original problem. This approach is 
outlined in Figure 7. 
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OPTIMIZATION METHOD 

MULTIPLE CONSTRAINT 
FORMULATION 


W: ■ WEIGHT OF i th SHIELD REGION 
D k ■ RADIATION LEVEL AT k th CRITICAL LOCATION 
C k - CONSTRAINT ON RADIATION LEVEL AT k th LOCATION 
A k - LAGRANGE MULTIPLIER 

tj ■ CHARAC TERISTIC THICKNESS OF i th SHIELD REGION 

1. DEFINE (INEQUALITY) CONSTRAINTS 
G k «D k <n-c k 

2. LAGRANGE FORMULATION IS. 

V(T,X) * W(Tl + XjGj - + * * ' + N( G k 

9 v av 

dtj ’ 0 ; ax k * °k 

3. FIND THE SET (I X > WHICH IS A SADDLE POINT OF V 

MIN MAX V (f£) 

tao X>0 

USING ARROW AND HURWICZ ITERATIVE PROCEDURE 

av av 
' IgRAcTvI 2 9t i 


av av 
Igrad vl 2 dX k 
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Optimality of Preferred Method 


OPTIMIZATION METHOD 
SNGLE CONSTRAINT 
FORMULATION 


Wj ■ WEIGHT OF SHIELD REGION 
1^ ■ RADIATION LEWI AT CRITICAL LOCATION 
C* • CONSTRAINT ON RAD IATION LEVEL AT k m LOCATION 
A - LAGRANGE MULTIPLIER 
n- EXPONENT SUFFICIENTLY LARGE TO EMPHASIZE 
DOMINANT CONSTRAINTS 

t f • CHARACTERISTIC THICKNESS OF i^ 1 SHIELD REGION 
1. COMBINE CONSTRAINTS INTO SINGLE RELATIONSHIP 



2. LAGRANGE FORMULAE ION IS 

V • Will + AUffi ; -&L ■ 0 
dt| 


3. SINCE THERE IS ONLY ONE LAGRANGE MULTIPLIER I A), 
FIND THE SET (T) SUCH THAT 



THEREFORE, REPETITIVELY MODIFY t m IN THE 
SET m WHERE m IS SUCH THAT 


L - WAX 


no aU/at m . MIN dUla *i 
' sm \ 0R aw/e^ 8W/at { 


FIGURE 7 


Experience with these methods on a variety of 
shield optimization problems led to a prefarence 
for the latter, single- constraint approach. One 
significant advantage was due to altering only one 
thickness value at each step of the iterative pro- 
cess, because this approach then reduced the 
amount of data that had to be updated prior to the 
next step. Hence, the single constraint approach 
was generally more efficient. It is likely that 
this advantage would be attainable in the multiple 
constraint approach if only the few largest Atj. 
were implemented at each step. 

The major objection to the conventional 
multiple constraint formulation lay in the fact 
that, for shielding problems, the measure of 
importance of a constraint is not fairly repre- 
sented by its linear distance from the criterion 
(as indicated in Item (1) of Figure 2). The 
exponential nature of the constraint functions 
results in the importance of each constraint being 

much more accurately expressed by the ratio of 
current value to criterion (as is utilized in the 
single constraint formulation). This failure to 
incorporate the exponential character of the 
function leads to severe difficulties in conver- 
ging on an optimum solution. The rate of con- 
vergence is necessarily dependent on the step 
size (i.e., AV) at each iteration and this is 
limited by the range of (At^, AX^) over which the 
partial derivatives can be applied. For the linear 
formulation indicated in the figure, this range is 
quite small. While some alteration in the expres- 
sion of the constraining relationships would 
presumably alleviate this difficulty, it was not 
pursued. 


The capacity of the single constraint method 
to produce a result which not only meets the 
criteria at low total weight, but is truly optimum, 
has been demonstrated numerically. The criteria 
for optimality are indicated for the conventional 
multiple constraint formulation in Figure 6, i.e.. 


3V 3W 
3t i“ 3t i 


3D, 


+ X 


k at. 


3D, 


i at. 


2 3t, 


0 ; 


i - 1, 2^.. I 


These equations can be evaluated using the set of 
thicknesses , t^, provided by the single constraint 
method directly and using a set of multipliers, 

Xfc, constructed from other data determined in the 


analysis : 


v n-1 


nXi / 

M . 3w/3t i 

\ 

V 

C J ’ 1 ’ 


There obviously are a total of I values of 
each X k , these values agreeing with one another 
to the extent that the values of ^ are in agree- 
ment. Since the essence of the single-constraint 
formulation is to align the values X^, the several 
values for each X^ can be made to agree quite well 
by aligning the X t within a very small difference. 
This is accomplished by continually reducing the 
increment (or decrement) to each variable once all 
constraints have been satisfied, iterating in the 
neighborhood of the optimum solution until satis- 
factory convergence is attained. This procedure, 
however, can be extremely time consuming and la 
of almost no benefit in terms of reducing the 
total shield weight from the value determined 
with a relatively coarse optimization criterion. 

Typically, the convergence of thickness values to 
the extent required to prove optimality doubles 
or triples the time required to determine a shield 
condition to within one percent of that optimum 
weight. 


SAMPLE OPTIMIZATION ANALYSIS 

c * C8UltS from a slm P lif ied application of the 
& WORD program are presented in Figure 8 to illus- 
trate some of the points made in the preceding 
discussion. This sample problem involves optimi- 
zation of the shield system Illustrated in Figure 

5 • 


Rather than modify the multiple constraint 
formulation to achieve a more effective technique 
for the kinda of shield optimization analyses of 
interest, the activity centered instead on exploi- 
tation of the single constraint approach. This 
latter technique, having fewer variables, was 
considerably simpler to apply and, where compar- 
isons were made of results obtained with both 
techniques, provided results that satisfied tests 
for optimality. 


The problem Included constraints on acute 
dose from solar cosmic rays (SCR) and on protracted 
dose from gecmagnetlcally trapped radiations and 
their secondaries, and from neutrons and gamma 
rays from an on-board isotopic power system. The 
total number of dose constraints was six, these 
being separate acute and protracted dose criteria 
for each of three critical organs: lens of the 

«ye. skin, and blood forming organs. 
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There were three man-model locations , repre- 
senting crew stations in a command area, an 
experimental area and a living area. The command 
area was specified to be within a biowell for SCR 
protection, and all SCR exposure was assumed to 
be taken at this location. Two of the man-model 
locations were on centerline and, with the vehicle 
being reasonably symmetrical cylindrically, a 
single dose point in the man-model was used to 
represent the location of each of the three organs. 
The man-model position in the biowell, however, 
was sufficiently off axis that some recognition 
of the resultant asymmetry was required. This was 
accomplished by using two dose points to represent 
each organ, one at the appropriate depth below the 
surface of the skin facing the outer wall and the 
other 180° opposite. The time-weights associated 
with the biowell locations were halved to compen- 
sate for the extra dose point. 

The value of making even token recognition of 
system asymmetry has been proven in a number of 
shield optimization problems. When dose can be 
delivered asymmetrically, as at the biowell location, 
it is particularly important that the dose point (s) 
be representative; otherwise. Introduction of an 
appreciable bias in the discrete and discontinuous 
representation of critical organ location can lead 
to false concentrations of shield mass. Ideally, 
the dose would be evaluated by sampling from a 
continuous timeline of astronaut position within 
the vehicle, sampling also from the solid angle 
at each position along the timeline; the effects 
of astronaut orientation at each position might 
also be diminished, since these are not necessarily 
real. This approach to dose integration has been 
incorporated in a recent version of the SIGMA 
program, as reported in Reference 6, but is not 
yet incorporated in SWORD. 

There were eight variables, three of which 
defined a cylindrically symmetric shield system 
enclosing all Inhabited areas, and five of which 
defined the biowell. 

Dose transmission data were furnished for a 
synchronous orbit mission of several months 
duration, these then being trapped electron and 
secondary bremsst rah lung dose, and SCR dose. The 
geometry, materials, and radiation source charac- 
teristics for the on-board isotopic power system 
were also defined and the direct and scattered 
doses computed therefrom. Only the direct dose 
proved to be of any significance and it had 
relatively little influence on the optimum shield 
determination, because this was dominated by SCR 
considerations; i.e., the acute dose criteria 
were the most stringent for the particular mission. 

Mission-integrated critical organ doses were 
calculated using a relatively coarse integration 
grid, namely four intervals in m (cosine of polar 
angle) and six intervals in 0 (azimuthal angle) 
or a total of 21 rays for each of the twelve dose 
points. The 252 rays traced made a total of 382 
shield crossings because of the nesting of the 
biowell shield within the main shield. The 
computer time required to perform the initial 
dose evaluation, in preparation for the shield 
optimization computation, was 21 seconds on the 
CDC 6500 computer. (Much more detailed integra- 
tion grids can be used; problems involving up to 
3000 rays have been solved.) 


The optimization history for this problem is 
shown in Figure 8. A total of 106 iterations on 
weight were accomplished in 40 seconds, each 
iteration initially involving weight increments 
of up to 250 Kg, decreasing to 15 Kg as the 
solution converged. It can be seen that two of 
the constraints dominated the optimization, the 
domination being particularly evident when the 
ratios of dose to constraint value decreased to 
approximately 0.9 after iteration 80, and when 
these ratios were then raised to an exponent of 
100. As noted on the figure, essentially the 
same result was obtained when the exponent was 4 
rather than 100; there is no particular reason 
not to use a large value of the exponent, however, 
other than to avoid overflow when the ratios are 
large, and SWORD internally increases the value 
to 100 as rapidly as the values of the ratios 
permit. 


TYPICAL SHIELD OPTIMIZATION PROBLEM 1<? 559 

(ITERATION HISTORY) 


• SIX DOSE CRITERIA 

• THREE TIME-WEIGHTED MAN MODEL 
LOCATIONS 


• EIGHT SHIELD THICKNESS VARIABLES, 
FIVE DEFINING BIOWELL 

•ONBOARD ISOTOPIC POWER SYSTEM 
•DOSE RATIO EXPONENT = 100 

•TWO CRITERIA DOMINATE; IDENTICAL 
RESULT WHEN FOUR IN EQUAL I TV 
CONSTRAINTS DELETED 

•DOSE EVALUATED WITH 252 RAYS 
IN 21 SEC (CDC 65001 


•OPTIMUM CONDITION DETERMINED 
WITH WEIGHT INCREMENT OF 
250 kg IN 40 SEC 


• WEIGHT INCREASED 5 PERCENT 
WHEN DOSE RATIO EXPONENT = 4 



FIGURES 


The optimization of the solution was checked 
at the 106th iteration by calculating the eight 
possible values (i.e., eight variables) of each 
of the six Lagrange multipliers (i.e., six dose 
criteria) from the several values of the multi- 
plier determined in the single constraint solu- 
tion. The set of multipliers obtained from each 
of the I equations was used to evaluate all I 
equations. Actually, each equation was divided 
through by A* so that the magnitude of the 
remainder was referenced to the magnitude of the 
numbers involved, i.e.. 
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The iteration was then continued for another 41 
steps, progressively decreasing the weight incre- 
ment to 1 Kg, with these results then being 
compared to the optimality as evaluated after 106 
steps. While the total shield weight decreased 
only 0.14%, the small changes in values of each 
of the variables led to a dramatic change in 
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apparent optimality: the A lk for the two .ignif- 

Icant constraints agreed within 0.4Z and the 
maximum value of R* for any set (A lk ) was 0.001. 
This experience is typical and serves to illustrate 
that the shield mass converges much more rapidly 
than the variables or the multipliers, and that 
the convergence criterion for iteration cutoff 
should be based on mass (as it is) rather than on 
some optimization test. 


CONCLUSIONS 

The work reported herein was initiated when 
a need arose to analyze the dependence of space 
vehicle shield weight requirements on various 
system parameters. The space vehicle in question 
was to fly long duration missions in synchronous 
orbit and hence a considerable mass of supplementary 
shielding appeared necessary. In order for the 
shield parametric data to be sufficiently self- 
consistent that the dependence on system parameters 
be meaningful, it was evident that a systematic 
procedure was required to produce it. The optimi- 
zation technique, and the implementation thereof, 
as described, proved to be extremely effective, 
both in terms of the scope of the problems which 
could be analyzed and optimized, and in terms of 
the significance of the results. Results from 
some production applications of SWORD were reported 
in References 15 and 16. 

The applicability of the optimization approach 
embodied in the single constraining function is 
much broader than that reviewed here. It is 
certainly worthy of consideration for, and has in 
fact been applied to, other shielding optimization 
problems. It could also be effectively applied 
in other, non-shielding efforts, possibly with 
some modification of the constraining functions 
to best approximate the functional relationships 
of constrained quantities to their specific 
constraint values. 
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